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PREFACE 

This  thesis  is  a  description  of  WEED,  an  algorithm  r 
simplifying  systems  of  equations  symbolically  and  classifying  variables 

r  further  numerical        >.   It  J  iphics- 

Oriented  Simulation  and  I-bdeling  System  currently  being  developed  at  the 
University  of  Illinois. 

Through  a  consistent  method  of  elimination,  WEED  produces  a 
condensed  system  of  equations  and  an  Equivalence  Table  for  variables, 
which  allow  for  a  more  efficient  compilation  and  numerical  integration 
than  would  otherwise  be  possible. 
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1.   THE  MODI.      YSTEM 

The  Geno  Lon  and  Modeling  System  currently  being 

developed  at  the  University  of  Illinois  has  as  its  primary  objective 
"to  pr      a  tool  which  he  can  (graphically) 

construct  models  of  a  proposed  system  and  analyze  them  for  static, 
dynamic  and  oscillatory  behavior.   The  user  will  have  the  capability 
to  define  elements  which  are  models  of  individual  items  of  his  system 
(e.g.  diodes,  motors,  water  pumps)  and  specify  their  behavior  by  means 
of  differential  and  algebraic  equations.   The  elements  will  have 
terminals  (e.g.  wires,  mechanical  links,  pipes)  through  which  certain 
variables  may  be  transmitted  (e.g.  voltage,  current,  torque,  displacement, 
pressure,  flow).   The  terminals  from  several  elements  can  be  connected 
together  at  nodes  (connection  points)  to  form  a  network. "   A  network 
may  be  used  as  an  element  of  a  larger  network,  and  hence  in  a  recursive 
fashion  fairly  complex  networks  can  be  constructed  quite  easily.   The 
overall  network  is  then  the  model  of  the  proposed  system,  and  is  equivalent 
to  the  system  of  differential  and  algebraic  equations  derived  from  each 
element  and  node  of  the  model.   Thus  the  model  can  be  completely  analyzed 
by  solving  this  equivalent  system  of  equations. 

The  user  interacts  with  the  Modeling  System  through  a  drawing 
program  which  builds  data  structures  representing  each  element  and 
network.   From  these,  Item  and  Global  Analyses  create  new  data  structures 
which  no  longer  contain  graphical  information.   The  equations  are  parsed 
into  tree  structures,  assigned  values  are  substituted  in  the  equations, 
and  all  external  names  are  replaced  with  new  internal  pointers.   WEED 
then  simplifies  the  equations,  classifies  them  according  to  the  types  of 
operations  and  variables  they  contain,  and  weeds  out  those  which  are  no 


2 

longer  needed.  Variables  are  organized  in  an  Equivalence  Table  according 
to  their  original  definition,  and  placed  into  different  sets  depending  on 
how  they  are  used  in  the  equations.  New  internal  names  are  then  assigned 
to  all  the  variables,  and  *„he  equations  are  updated  and  restructured  into 
a  compact  format.  The  final  stage  then  compiles  the  condensed  system  of 
equations,   and  solves    for  transient,   steady-state  or  oscillatory   response. 
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2.      PRELIMINARIES   AND  OVERVIEW   OF  TOE  ALGORITHM 
Variables    contained  in   the  system  of  equations    are  not   all  of 
the   same  nature.      Some   refer  to  terminals   in   the  network,    and  are 
identified  as    I -variables    (e.g.    electrical  current).      Others   are 
associated  with  nodes   and  are   classed  as   E-variables    (e.g.    voltage). 
iOcai  variables    (also  called  internal  variables)    are   those  which   are 
specific   to   an   instance   of  element   in   the  network.      Global  variables, 
on   the   other  hand,    are    common  to   groups   of  elements    (e.g.    ambient 
temperature).      The   special   global  variable   TIME  is  permanently  assigned. 
Finally,   parameters   are   those  variables   associated  with  a  given   type   of 
element,  whose  value   is    specified  with  each   instance   of  that   element   in 
the  network    (e.g.    resistance   of  a  resistor). 

In  addition  to  keeping  track  of  the  above  categories,  WEED 
places  all  variables  into  four  mutually  exclusive  sets  ,  according  to 
their  usage   in  the   equations : 

a.  Set  S — all  global  variables   and  variables    defined  by  S   - 
equations    (of  the    form  V  =    function    (globals   and 
constants   only)). 

b.  Set   S„ — variables    defined  by  S   -equations    (of  the   form 
V  =   function    (globals,    constants,    and  time)),    and  the 
special  variable   TIME. 

c.  Set   L — variables   involved  only   in   linear  expressions. 

d.  Set  M — variables   involved  in  non-linear  or   differential 
expressions.      As   an   example,   given   the  equations 


EQ  =   3*E1  *  G0*(E2+E3)    -   1/E^   +    (E5+2)»T  -  Eg«E     +££ 

h  "  3  -  V2 
lx  =  I0*sin(   *T) 

h  -   <V3)/<IlV 

where  the  E.    are   E-variables,   the   I.    are   I -variables ,    the  G.    are   globals 
111 

and  T   is   tine,   E     through  E     would  belong   in   set   L;   E,    through  Eo   in   set 
M;   I      in   set  S    ;    and  I     and  I      in  set  S    .      Note  that   since   global  values 
are   fixed  throughout  any  given  execution   of  the  Modeling  System,   global 
variables   are   treated  as    constants   insofar  as   variable   and  expression 
classification  is   concerned. 

An  additional   set   DY  of  new  variables   is   created  for  the 
output   of  general  equations,   and  is   discussed  in   section  h. 

A  number  of  passes   are  made  through  the  system  of  equations, 
each  time   advantage  being  taken  of  simplifications   made   in  previous 
passes.      As   each  equation   is  processed,   all  expressions  which   are   only 
functions   of   constants   are   evaluated  immediately  and  replaced  by  the 
result.      Expressions   of  the   form 

E+0,    0+E,   E*0,    0*/E,   E*/l,    1*E,   E**0,   E**l, 

1**E,  and  —  (S.,  or  S  -expression) 
at   l     d 

are  simplified  as  they  are  detected,  and  unary  minuses  are  eliminated 

wherever  possible.   Equations  which  become  trivial  relations,  such  as 

variable  =  constant 
or 

variable  =  +_  variable 
are  eliminated  by  setting  appropriate  pointers  in  the  Equivalence  Table 
(described  in  section  U),  and  the  remaining  equations  are  placed 


into  three  sets:   S  -equations,  S  -equations,  and  general  equations 
(those  which  are  neither  S  nor  S  ). 

The  reason  for  arranging  all  the  equations  and  variables 
into  different  sets  is  that  each  set  can  receive  special  treatment  in 
the  numerical  analysis  phase,  thereby  increasing  performance  and 
reducing  cost.   For  instance,  variables  in  set  S  need  only  be  evaluated 
once  per  simulation;  those  in  set  Sp  once  per  time  step;  and  derivatives 
of  variables  in  set  L  need  not  be  evaluated  at  all. 

The  equations  are  reprocessed  until  there  are  no  more  changes 
in  the  Equivalence  Table  (up  to  a  maximum  of  three  passes).   All 
variables  are  then  renumbered  so  that  the  members  of  each  set  are 

contiguous  in  memory — that  is,  the  numbers  0  to  n   -1  refer  to  the  n 

Sl  Sl 

members   of  set  S    ;   n        to  n        +  n      -1  to  the  n        members   of  set  S    ; 
1        S1  S1  b2  b2  2 

n       +  n        to  n        +  n        +  n   -1  to   the  n     members   of  set   L;   and  n        +  n 
Sl  S2  Sl  S2  L  L  Sl  b2 

+  n     ton        +n        +nT+  n,  -1  to  the  n„  members   of  set  M.      The  new 

L     S     S     L    M  M 

names  are  substituted  in  the  equations,  and  these  in  turn  are  compressed 
and  output  in  their  three  different  sets  to  the  numerical  integration 
routines.   A  concise  flow  chart  of  the  WEED  algorithm  is  given  in 
Appendix  A. 


3.   INPUT  SPECIFICATIONS 
The  input  to  WEED  consists  of  a  set  of  Input  Tables  and 
a  single  list  of  equations,  as  shown  in  Figure  1. 


fullvord- 


-   ptr.    const,    table 


L'    all   Input 


ptr.    equation   list 


#  E-variables 


#   I-variables 
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location  In 
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name  table 
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Equation   1 
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Equation  2 
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displacements  from 
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doubleword ■ 


Figure  1.   Input  to  WEED 


7 
The  Name   Tables    contain   the  actual  eight-character  variable 
names   provided  by  the  user.      In   the  equations,   however,    these  variables 
are   represented  by  the   displacement   of  their  name    from  the  beginning 
of  their  respective  Name   Tables.        The  same   holds   true    for   constants, 
where  the   Constant   Table   contains  the  double-precision   floating-point 
constant  values. 

The  equations   are   in  tree-structured  form,    and  consist  of  a 
series   of  "nodes"   of  information   connected  by  pointers,   as   illustrated 
in  Figure  2.      The   first  node   of  each  equation   consists   of  a  pointer  to 
the  next   equation   in   the   list,    as   a  displacement   from  the   current 

equation.      The   remaining  nodes   are   clusters   of  halfwords    (usually  two 

2 
or  three)      contiguous    in  memory.      The   first  HW  of  each  node   is   an  op-code 

which  indicates  the  type  of  node.      If  the  node   represents   an   operation, 

3 

the  remaining  HW- fields  contain  pointers  to  its  respective  operands. 

If  the  node  is  a  variable  or  constant,  the  second  HW  contains  the 
internal  name  associated  with  it.   Table  1  describes  the  different  op- 
codes and  associated  operands. 


E-  and  I-type  variables  do  not  have  alphameric  names  associated  with 
them,  and  are  numbered  sequentially,  starting  from  0. 

2 

The    'user-specified  function'    node   is   a  variable-length  node, 

including   fields    for   function   type   (second  HW)    and  number  of 

pointers    (third  HW). 

3 

A  'function'  node  contains  a  field  specifying  the  type  of  function 

(second  HW). 
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Figure  2.   Equation  Format 
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Internal   Codes    for  Equation  Nodes 


10 
•  •   OUTPUT  REQUIREMENTS 
The  output  consists  of  a  Table  Block  and  an  Equation  Block, 
iepicted  Li    ;  ;re  3.   The  Table  Block  consists  of  the  Constant 
Tab.  ig  possl-      litions  create:  -    le  WEED  step),  the  input 

•  use  i.      3-referencing  and  diagnostics),  and  the 
.  le  (E'v,   .    ich  contains  all  the  variables  ordered 
according  to  network  type,  with  their  respective  new  internal  names 
(variables  found  to  be  equivalent  will,  of  course,  have  the  same  assigned 
name).   The  Table  Block  also  contains  other  necessary  data,  such  as 
rs  to  the  various  tables  and  to  the  Equation  Block. 

Upon  entry  to  WEED,  the  variables  are  arranged  in  the  EQV  in 
order  of  ascending  numeric  names  for  each  field  in  the  table,  as  shown 
in  Figure  k.      As  the  equations  are  processed,  different  types  are 
assigned  to  the  variables,  and  appropriate  pointers  are  set  for  those 
which  are  found  to  be  equivalent  to  other  variables  or  to  constants. 
During  internal  name  reassignment,  all  variables  are  marked  'assigned' 
(except  for  those  equivalent  to  constants),  so  that  upon  output,  only 
three  types  will  be  present  in  the  EQV:   assigned  (+V),  assigned  (-V), 
and  constants. 
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TYPICAL  CONFIGURATION 
TYPE     VALUE 


EQV  TYPE 

REPRESENTATION 
1st  HW    2nd  HW 
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X'OO' 
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in 
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CI 
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const,  table 

+ Variable 

Fl 

ptr.  to  var 
in  EQV 

Displ. 

in 
Table 

-Variable 

FF 

ptr.  to  var 
in  EQV 

GL0BALS        ( 


nsra 

VARIABLES 


parameters/ 


- 


I-VARS 


<«  .alfword-*.-«Halfword-». 

Figure   k.      Equivalence   Table 
The  Equation  Block   consists   of  three   sub-blocks    corresponding 
to  the  three   different   sets   of  equations.      The  equations   take   a  form 
similar  to  the   input   equations,  with  two   differences: 

1.      Operation  nodes   incorporate   operands  which   are  variables, 
their  pointer  fields  being  replaced  by  the   corresponding 
new  variable  names.      The  remaining  pointers   to  other  nodes 
take   on   a  negative  value,   in  order  to   differentiate    from 
variable  names. 


13 
General  equations  are  put  in  the  form  0  =  expression, 
and  the  0  is  replaced  by  a  new  variable  of  the  set  DY, 

where  the  DY's  are  numbered  starting  at  nc   +  n^  +  n   +  n 

Sl    S2    L 

. -gure  5  illustrates  the  modifications  performed  on  the 
example  of  Figure  2. 
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Figure   5-      Output  Equations 
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>.       EQUATION  CLASSIFICATION 

For  this  and  the  following  sections,  the  flow  charts  provided 
in  the  Appendices  should  be  useful  as  an  aid  to  understanding  the 
..m. 

Upon  entry,    3  goes  through  an  initialization  phase,  in  which 
s  an  area  for  the  output  Table  Block,  copies  the  Name  Tables 
and  Constant  Table  from  the  input  area,  sets  pointers  to  the  different 
tables,  and  zeroes  out  the  Equivalence  Table.   It  is  then  ready  for  the 
Equation  Classification  stage. 

WEED  takes  each  equation  in  turn,  replaces  the  '='  in  the  first 
node  by  a  '-',  and  calls  upon  the  PROCESS  Routine  (described  in  section  6  ) 
to  determine  the  types  of  expressions  on  the  left  and  righthand  sides 
of  the  equation.  PROCESS  travels  recursively  down  the  tree,  performing 
simplifications  as  it  goes  along,  and  returns  two  pairs  of  numbers 
representing  the  left  and  right  expression  types.   Every  expression  has 
two  numbers  associated  with  it:   one  (referred  to  in  the  program  as  EXPTP, 
LX  or  RX)  which  gives  an  indication  of  the  structure  and  complexity  of 
the  expression;  and  another  (referred  to  as  S1S2EXP,  S1S2LX,  or  S1S2RX) 
which  indicates  whether  the  expression  is  type  S  ,  S  ,  or  neither. 
Table  2  lists  the  possible  types  and  respective  codes. 
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(EXPTP,  LX,  RX)  (S1S2EXP,  S1S2LX,  S1S2RX) 

0  constant  '0*     (S1S2EXP=1)  0  Neither  S  nor  S  -expression 

1  other  constant           =l)  1        ression  (single  consts.  included) 
single  variable  2  S  -expression 

3  constant  +_  variable 
h  variable  +  variable 
5   general  expression 


Table  2.   Expression  Types 

Once  WEED  receives  the  left  and  right  expression  types  from 
PROCESS,  it  restores  the  '='  and  tries  to  put  the  equation  in  a  standard 
format,  according  to  the  following  priorities: 

1.  S  or  S  -expressions  on  the  right  side.   (if  the  left 
side  is  also  an  S  or  S  -expression  an  error  is  indicated, 
the  equation  is  eliminated  and  the  next  equation  is  fetched. ) 

2.  Single  variable  on  the  left  side. 

3.  Doublet  (C+V  or  +V  +V)  on  the  left  side. 

After  the  appropriate  switches  are  made,  if  a  doublet  (or  any 
expression  of  the  form  E+E)  is  found  on  the  left,  it  is  split  up,  the 
left  branch  of  the  doublet  being  transposed  to  the  right  side  of  the 
equation.    At  this  stage,  the  equation  will  necessarily  be  in  one  of 
the  following  forms : 


Any  constants  in  a  doublet  are  always  kept  on  the  left.   Hence  in  a 
split-up,  the  constant  always  gets  moved  to  the  right  side  of  the  equation, 


Lable  =  +      ole 
variable  =  expresc  i 
- 
We  scribe  th  .  each  case. 

5-1  Variable  =  Const 

Set  variable  type  to  'constant'  (X'Cl')  in  the  EQV.   Place 
a  pointer  to  the  constant  (as  a  displacement  from  the  beginning  of  the 
Constant  Table)  in  the  EQV  value  field  of  the  variable.   Eliminate  the 
equation,  and  set         to  1,  indicating  a  change  made  to  the  EQV. 


p. 2   Variable  =  Variable 

The  two  variable  types  determine  the  action  to  be  taken, 
according  to  the  following  table: 

Variable^ 


Variable 


- 


L 

C 
+  7 


Sl   S2   L 


+V   -V 


E 

E 

£ 

-*• 

E 
1 

_3k 

_^ 

E 

E 

^- 

^- 

i 

1 

E   1 

1 

T— 

OK 

OK 

OK 

OK 

ok| 

OK 

OK 

OK 

OK 

— 

OK 

OK 

OK 

OK 

E 

~1- 

_i 

E   1 

I 

— 

"> 

^ 

~^- 

N 

^ 

^ 

. 

'% 

i- 

_i 

">  j 

_ 

♦  ' 

IT 
E:   Error.   Relay  message  to  the  user  and  delete  the  equation. 
4^:   Switch  the  variables,  so  that  variable  becomes  variable 

and  vice  versa.   Refer  to  the  table. 

2 
^  :   Follow  the  pointer  (or  chain  of  pc     ■  ;)   beginning  at 

■iable,  .   The  last  variable  in  the  chain  becomes  the 

new  variable  .   Return  to  the  table. 

OK:   Set  the  type  of  variable  to  +V  or  -V  according  to  the 

equation  which  defines  it.   If  variable,,  points  to  another 

d. 

2 
variable,   follow  the  pointer  chain.        Place   a  pointer  to 

3 
the  new  variable      in   the  value    field  of  variable    . 

Eliminate   the   equation   and  set  EQVFLAG   to  1. 


5 . 3      Variable  =   Expression 

If  the  expression   is   neither  a  S-nor  S   -expression,   go  to 
section  5.*+    (Expression  =  Expression).      Set   the  variable   type  to  S 
or  S     accordingly.      If  the  equation  has   already  been  placed  in  the 
proper  equation   set,    proceed  to   the  next   equation.      Otherwise,   mark 

:   equation  by  placing  a  1  or  2   in  the   first  byte   of  its   link   field, 
and  set  EQVFLAG  to  1. 


.  order  to  avoid  following  long  chains  in  future  searches  through 
the  EQV,  a  direct  pointer  is  set  from  the  variable  being  traced  to 
the  last  variable   in   the   chain. 

- 

-  .  the  case  where  variable^  points  to  a  cor..  *  n1  ,  variable   is  made 

to  point  to  the  same  constant  (or  its  complement)  and  its  type  is 
set  to  'constant.' 
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j . k     expression  =  Expression 

If  the  equation  is  in  any  of  the  following  forms,  simplify 
it,  retest  it  for  further  simplifications  and  reprocess  it: 


VE2  " 

0 

-> 

El 

— 

0 

EXC  = 

0 

-> 

El 

= 

0 

for  C  4   0 

CE1  * 

0 

->■ 

El 

= 

0 

for  C  i   0 

El  V 

0 

-> 

El 

= 

0 

EliE2  = 

°3 

-> 

El 

= 

"E2 

El/E2  = 

E3 

->■ 

El 

= 

E2E3 

E3    * 

El/E2 

->■ 

El 

= 

E2E3 

El/E2  "  E3/E^  "  El  \-   E2  E3 

If  from  the  first  test  it  is  found  that  the  equation  fits  none 
of  the  above  categories,  place  it  in  the  set  of  general  equations  (by 
putting  a  3  in  the  first  byte  of  its  link  field)  and  proceed  to  the 
next  equation. 
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6.   THE  PROCESS  ROUTINE 

Given  a  pointer  to  the  top  node  of  an  expression,  PROCESS  will 
simplify  the  expression  and  determine  its  type.   The  types  of  lower  level 
expressions  (i.e.  operand  expressions)  are  used  in  determining  the  type 
of  the  current  expression.   Thus  pairs  of  numbers  (EXPTP ,S1S2EXP) 
representing  each  lower  level  expression  are  transmitted  recursively, 
from  the  bottom  up. 

The  top  node  of  each  expression  determines  the  subsequent 
steps  to  be  taken,  as  follows: 

6  .1    ..  -.  ;.-. 

Determine  whether  the  constant  is  zero  or  non-zero,  and  set 
EXPTP  to  0  or  1  accordingly.   Set  S1S2EXP  to  1.   Return. 

6.2   Global  Variable  or  TIME 

Set  EQV  type  to  S  if  node  is  a  global,  or  to  S  if  it  is 
TIME.   Go  to  5.3. 


6 . 3  Other  Variable 

If  unused  so  far,  make  the  variable  type  L.   If  type  L  or  M, 
set  S1S2EXP  to  0.   If  type  S   or  S   set  S1S2EXP  to  1  or  2,  respectively. 
If  type  +V  or  -V,  follow  pointer  chain  and  repeat  5.3.   If  type 
'constant,'  replace  node  in  tree  with  a  corresponding  constant  node,  and 
go  to  5.1.   Set  EXPTP  to  2  (single  variable)  and  return. 
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Binary  Operation 


Process  the  left  and  right  expressions,  remembering  their 
respective  types.  S1S2EXP  for  the  current  expression  is  determined 
from  the  following  table : 

S1S2RX  (right) 


S1S2EXP 

0 

1 

2 

0 

0 

0 

0 

1 

0 

1 

2 

2 

0 

2 

2 

S1S2LX 
(left) 


Simplify  unary  minuses  according  to  the  rules  of  algebra,  placing  any 
remaining  unary  minus  above  the  current  node  (in  most  cases,  unary  minuses 
can  be  successively  bubbled  up  in  this  manner,  so  that  upon  reaching 
the  top  of  the  equation  at  most  one  unary  minus  remains). 

If  the  operation  is  addition  or  subtraction,  eliminate  any 
zeros  on  either  side  of  the  expression.   Add  or  subtract  constants 
appearing  in  constant  nodes  or  doublet  expressions  (C+V),  keeping  the 
resultant  constant  on  the  left  branch  (either  of  the  current  node  or 
of  the  left  doublet).   Determine  EXPTP  from  the  left  and  right  expression 
types  according  to  the  table  below,  and  return. 


1 


If  the  top  node  of  the  equation  is  being  processed  (with  '  =  '  replaced 
by  '-♦)  zeros  are  preserved. 
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RX 


DC 


EXPTP 
(+/-) 

0 

1 

2 

3 

1* 

5 

0 

1 

2 

3 

It 

5 

1 

1 

1 

3 

3 

5 

5 

2 

2 

"3 

i 

*5 

5 

3 

3 

3 

5 

5 

5 

5 

U 

it 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

5 

If  the  operation   is  multiplication  or  division,   eliminate   any 

l's   or   -l's   on  either   side   of  the  expression.      If  either  side   is   0,   make 

2 
the   current  node   a  0.        If  both   sides  are   constants,   evaluate   the 

expression.      From  S1S2LX  and  S1S2RX   determine  whether  the  expression   is 

non-linear  under  multiplication   or   division,    as    follows: 

S1S2RX 


CALL 
NONLIN 

0    0    2 

0 

*,/       •■/ 

1 

/ 

2 

*,/ 

S1S2LX 


If  it   is  non-linear,    call  NONLIN   (a  routine  which  sets   all   L-type   variables 
in  the   subtree   to   M-type   and  marks   the   top  node   of  the  expression   as 
non-linear  by  making  its   op-code  negative.      This  way,    if  a  higher  level 
expression   is   also   found  to  be  non-linear,   the  present   subtree  will  not 
have   to  be  rescanned) .      Finally,    determine  EXPTP   from  the  table  below, 
and  return. 


Division  by   zero  causes   an  error  message   and  automatic  elimination  of 
the  equation. 
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RX 


LX 


EXPTP 

0 

1 

o 

■  --  ■  i 
3 

k 

5 

(*,/) 

•1' 

or    '-1' 

other 

0 

0 

0 

0 

0 

0 

0 

0 

1 

'1* 
or   '1' 

u 

]_ 

1 

2 

3 

1+ 

5 

other 

0 

1 

1 

5 

5 

5 

5 

2 

0 

2 

5 

5 

5 

5 

5 

3 

0 

3 

5 

5 

5 

5 

5 

U 

Q 

It 

5 

5 

5 

5 

5 

5 

0 

5 

5 

5 

5 

5 

5 

For  the  case  of  a  power  operation,  if  both  operands  are 
constants,  perform  the  exponentiation  (C  **C  evaluated  as  e  ^    l)s 
and  set  EXPTP  to  0  or  1  depending  on  the  result.   Perform  the  following 
simplifications  if  possible: 

E**0  -»■  1   (EXPTP  <-  1) 

E**l  ■>  E   (EXPTP  *■   LX) 

0**E  +   0   (EXPTP  «-  0) 

1**E  ■+  1   (EXPTP  ♦■  1) 
If  no  simplification  is  possible,  call  NONLIN  and  set  EXPTP  to  5 
(general  expression).   Return. 

If  the  operation  is  max  or  min  evaluate  the  expression  if 
both  operands  are  constants.   Otherwise,  call  NONLIN  and  make 
EXPTP  =  5-   Return. 


6 .5   Unary  Operation 

Process  the  left  operand.   S1S2EXP  is  passed  from  operand  to 
operator  unchanged  (i.e.  31S2EXP  «-  S1S2LX) . 
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If  the  operation   is    differentiation,   move   any  unary  minus 
upwards.      If  the  operand  is   an  S   -expression,   replace  the   current 
expression  with  the   constant    '0'    and  set  EXPTP  to  0.      Otherwise, 
call  NONLIN   and  make  EXPTP  =   5.      Return. 

For  a  unary  minus   operation,    if  the  operand  starts  with  a 
unary  minus,    delete  both  unary  minuses.      If  the  operand  is   a 
constant,    complement   it    and  delete   the  unary  minus.      Copy  EXPTP   from 
the  operand   (i.e.    EXPTP  «-    LX)    and  return. 

If  the  operation   is   function ,   evaluate   it   if  the  operand  is 
a  constant,   replacing  the  expression  with  the   result,   and  setting 
EXPTP  to  0   or  1.      Otherwise,    call  NONLIN   and  set  EXPTP  to   5.      Return. 

6.6      User-Defined   Function 

Process   all  the   operand  subexpressions,    call  NONLIN,    set 
S1S2EXP  to  0   and  EXPTP  to  5-      Return. 
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1.       THE  FINAL  STAGES:   NAME  REASSIGNMENTS , 
EQUATION  RESTRUCTURING  AND  OUTPUT 

Having  processed  the  entire  system  of  equations,  WEED  counts 
the  numbers  of  independent  S  ,  S  ,  L  and  M  variables  in  the  EQV  and 

--.ializes  each  of  four  registers  with  the  first  numeric  name  to  be 
assigned  to  variables  in  each  of  the  four  sets.   As  the  EQV  is  scanned, 
each  time  a  new  variable  is  encountered  its  value  field  is  overwritten 
with  the  contents  of  its  associated  register,  its  type  set  to  ' assigned(+V) 
and  the  register  incremented  by  1.   Whenever  a  variable  type  +V  or  -V  is 
encountered,  the  new  name  of  its  equivalent  variable  is  copied  in  its 
value  field,  and  its  type  set  to  ' assigned(+V) '  or  'assigned(-V) ' , 
accordingly.   'Constant'  type  variables  remain  unchanged.   As  a  final 
result,  all  variables  will  have  been  numbered  sequentially,  with  all 
S  -variables  first,  S_-variables  next,  and  so  on. 

At  this  point,  the  equations  are  modified  so  that  the  new 
variable  names  are  contained  within  operation  nodes.   As  the  equations 
are  processed,  a  cumulative  byte  count  is  kept  for  each  of  the  three 
equation  sets ,  so  that  an  appropriate  amount  of  core  can  be  later 
requested  for  the  output  Equation  Block.   Each  tree  is  traversed 
recursively  by  means  of  the  TRAVEL  Routine  (MODFORM  mode).   As  each 
variable  is  detected,  it  is  incorporated  in  its  parent  node  by  over- 
writing the  respective  pointer  field.   An  exception  if  made  for  the 
case  where  the  variable  is  ' assigned(-V) ' ,  in  which  case  the  old  variable 
is  transformed  into  a  unary-  node  whose  operand  field  contains  the  new 
variable  name.   Constants  are  not  incorporated  in  parent  nodes,  and  all 
remaining  node-pointer  values  are  made  negative. 


25 
Once  all  the  equations  have  been  reformatted,  WEED  issues 
a  GETMAIN  for  the  output  Equation  Block  and  calls  upon  the  PUTOUT 
Routine  to  supervise  the  output  operation.   PUTOUT  looks  only  at  S  - 
equations  to  begin  with,  and  passes  them  one  by  one  to  TRAVEL  (TRANSFER 
mode),  which  recursively  moves  individual  nodes  to  consecutive  memory 
locations  in  the  output  area.   Since  the  equations  are  squeezed  into  a 
compact  form,  the  displacement  relationships  between  the  nodes  are 
altered,  and  hence  new  pointer  values  must  be  set  as  the  nodes  are 
moved.   New  link  pointers  to  the  next  equation  can  be  set  as  soon  as  each 
equation  transfer  is  completed. 

PUTOUT  repeats  this  process  for  S   and  General  Equations,  and 
returns  to  WEED  after  setting  pointers  to  the  three  sets  of  equations 
in  the  Table  Block.   Finally,  the  entire  input  area  is  FREEMAINed,  the 
compiler  is  called,  and  WE  ARE  DONE. 
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APPENDICES 
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I 


GIT  SPACE   TOR   OUTPUT 
TABLE   BLOCK,    COPY 
CONSTANT  I  NAME  TABLES, 
3ET  POINTERS  ,   ZERO  OUT 
EQUIVALENCE   TABLE 
CYCLES  -   0 


CYCLES  *    CYCLES   ♦   1 

EQVTLAG  -    0 

GET  INPUT  EQUATION   PTR 


GET  NEXT  EQUATION, 
CLASSIFY    IT   OR 

ELIMINATE    IT 


EQVFLAG  -    1 


SCAN  EQV  AND  COUNT 

#S1,   S2,    L1H  VARS 
ASSIGN  NEV 

INTERNAL  NAMES 


SUBSTITUTE   VARIABLE 


NAMES    IN  ALL  EQUATIONS 


k  MODIFY    FORMAT 


FOR  OUTPUT 


GET  SPACE   FOR  OUTPUT 
EQUATION   BLOCK.      MOVE 
31-EQHS  to  SI   BLOCK 
S2-EQNS  TO  S2   BLOCK 
GEN   EQNS   TO  GEN   BLOCK 
SET  PTRS    IN   TABLE  BLOCK 
FREEMAIN    INPUT  AREA 


SET  ALL  SI   i  S2   VARS 
(EXCEPT  GLOBALS   i  TIME) 
AND  ALL  M  VARIABLES  TO 
TYPE   L  IN   EQV 


Y 


N 


THE  WEED 
ALGORITHM 
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-jui   ai     •  an 

t  I  -■ 


-•/wnocau  V- 


nuitrow  to 
•  (un^  -  tin. 


UM,    •  1JML  -  IXP« 


VU  TYPE  »    CO»STMT 
3ET  PT»   TO 
COISTA»T  TABLE 


VAP  TYPE  -    S,    OR  S, 


ZQ)I    TYPE  -    S,    OB  S 


EQUATION 
CLASSIFICATION 


'-'   PVt 

•suiTTof  ;i 


CON    TYPE  *     J 

(CEXEHAL) 


DELETE 
OUOKIHATOP. 
Of*  EXPONENT 


MODTPY   E8UATI0H 
TO  HAKE 
OP   -   • 
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:■  mCBS 

LEFT   I    RI»T 
DPIMNOE 


crawcm 

SLS2IXP    FROM 
US2LX   4  S1S2KX 


■P  -    1 


Binary   Opirotlon 


SL32EXP  *    2 


S1S2EXPT  *   0 
EXPTP  *  5 
CALL  II  ON  LIN 


REPLACE    VARIABLE 
IN    TREE  BY 
CONSTANT 


VARIABLE   TYPE  *    L 


S1S2EXP  *    0 


_(J      RETURN 


Unary   Operation 


3IS2EXP*    1 


PROCESS 

LEPT 

OPERAND 


SLS2EXP  ♦   S1S2LX 


THE   PROCESS  ROUTINE 
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1 

Y 

r 

ADD  COKSTAXTS 
RESULT  TO   L  OB  (L.R) 
UPDATE    CONSTANT  TABU 
DEISTS   COHSTAIIT  ON 
P  OR   (L,Ri    RODE 

DETERMINE  EXPJ 
.''BOH   LX  AND  RX 


CALL  NONLIN 
EXPTP  *    5 


EVALUATE  EXPR 
UPDATE   OONST  TABLE 


REPLACE  EXPR 
WITH  RESULT 
EXPTP  »  0  OR  1 


REPLACE  EXPRESSION 
WITH   CONSTANT   '1' 
EXPTP  -   1 


DELETE   RIGHT 
NO  IE 
EXPTP  *    LX 


EVALUATE  EXPR  AND 
REPLACE  WITH    RESUUT 


EXPTP  -  0  OR   1 


CALL  NONLIN 
EXPTT  «    5 


PROCESS   ROUTINE 
(BINARY  OPERATIONS) 
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PROCESS  ROUTINE 
(MULT/DIV) 


MULTDIV 


I 


BUBBLE 
UNABY-'i 


' 

' 

OP  =   /?           ^ 

j£n 

CALL  NONLIH 

N^ 

■^S   -EXPR    ^V. 

^  ON    LEFT?  ^S 

*I 

EXPTP  -    5 

RETURN       14 
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:.   tast  'O1 


EVALUATE   FUNCTION 
REPLACE  EXPR  WITH 
RESULT.    UPDATE 
CONST.   TABLE 
EXPTP  *   0  OH   1 


BUBBLE 
UNARY- 
(IF  ANY) 


DELETE   BOTH 
UNARY-' s 


COMPLEMENT   CONST 
UPDATE   CONST 
TABLE    DELETE 
UNARY - 


:a^i  NONLIN 
'    -    5 


PROCESS   ROUTINE 
(UNARY  OPERATIONS) 
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